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Abstract. We report on results of two BeppoSAX observations of the Z source GX 17+2. In both cases the source is in the 
horizontal branch of the colour-intensity diagram. The persistent continuum can be fit by two-component models consisting of 
a blackbody plus a Comptonization spectrum. With one of these models, two solutions for the blackbody temperature of both 
the observed and seed photons for Comptonization are equally accepted by the data. In the first observation, when the source is 
on the left part of the horizontal branch, we observe a hard tail extending up to 120 keV, while in the second observation, when 
the source moves towards right in the same branch, the tail is no longer detected. The hard (> 30 keV) X-ray emission can 
be modeled either by a simple power-law with photon index F ~ 2.7, or assuming Comptonization of ~ 1 keV soft photons 
off a hybrid thermal plus non-thermal electron plasma. The spectral index of the non-thermal injected electrons is p ^ 1.7. 
The observation of hard X-ray emission only in the left part of the horizontal branch could be indicative of the presence of a 
threshold in the accretion rate above which the hard tail disappears. An emission line at 6.7 keV with equivalent width ~ 30 
eV is also found in both observations. We discuss these results and their physical implications. 
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1. Introduction 

Z sources form the brightest class of galactic Low Mass X-ray 
Binaries (LMXBs), with X-ray luminosities typically in the 
range ^ 0.5-1 times the Eddington luminosity (where ^Edd ~ 
1.5 M/Mq X 10'^* erg s~^). Their name arises from the pat- 
tern they trace in the colour-colour or hardness-intensity di- 
agram (CD and HID, respectively) on typical time scales of 
days. The three regions which characterize the Z-track, from 
the top, are the horizontal branch (HB), the normal branch (NB) 
and the flaring branch (FB). Multifrequency observations of 
ScoX-1 (Vrtilek et al. 1991a; Hertz et al. 1992; Augusteijn 
et al. 1992) and Cyg X-2 (van Paradaijs et al. 1990; Vrtilek 
et al 1990; Hasinger et al. 1990) have shown that the ac- 
cretion rate M increasing as the sources move in the direc- 
tion HB — > NB FB is consistent with UV observations 
and some QPO models. Several two-component models have 
been proposed to fit the low energy (^ 20 keV) X-ray spec- 
tra of Z sources, such as a simple or Comptonized blackbody 
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(bb) plus a multicolour disk blackbody (DBB, the "Eastern 
Model", Mitsuda et al. 1984, 1989), a BB plus an unsaturated 
Comptonization spectrum (F(£') oc E^^ exp {—E / Ec), the 
"Western Model", White et al. 1986, 1988) and, more recently, 
a DBB plus the Comptonization model worked out by Titarchuk 
(1994, COMPTT in XSPEC, Piraino et al. 2002; Di Salvo et al. 
2002). In recent years, transient hard X-ray tails extending up 
to 100 keV have been discovered from the classical Z sources 
(GX 5-1, Asai et al. 1994; Cyg X-2, Frontera et al. 1998; Di 
Salvo et al. 2002; GX 17+2, Di Salvo et al. 2000, hereafter 
DS2000; Sco X-1, D'Amico et al. 2001; GX 349+2, Di Salvo 
et al. 2001) and from the peculiar Z source Cir X-1 (laria et 
al. 2001). All these hard tails have been fit with power-laws 
(PL, F{E) oc E^^) with photon index F ranging from ~ — 1 
(Sco X-1, when the source was in the FB) to ~ 3.3 (Cir X-1). 
The time behaviour of the hard tail is complex and still not well 
understood. DS2000 observed a strong hard tail in GX 17+2 
when the source was on the left part of the HB; as the source 
moved on the right across the branch, an independent fit of the 
photon index and PL normalization was not possible, and fix- 
ing the photon index at 2.7 yielded an intensity of the PL that 
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lowered (chance probability of F-test ^ 10^'^) until it disap- 
peared at the HB-NB apex. In the case of Cyg X-2 (Di Salvo et 
al. 2002) the hard tail was observed when the source was in the 
HB, while in a later observation, when the source was in the 
NB, no evidence of hard X-ray emission was found. The hard 
tail of GX 5-1 (Asai et al. 1994) was detected when the source 
was in the NB, and decreased as the source moved to the FB. 
If M actually increases from HB to FB, these results could be 
indicative of the presence of a anti-correlation between hard 
X-ray emission and accretion rate M. However, the results of 
D' Amico et al. (2001) contrast with this scenario, as in Sco X-1 
the hard tail was observed in all branches without any apparent 
correlation with the source position on the CD. 

1.1. GX 17+2 

GX 17h-2 is one of strongest Galactic X-ray sources, and 
was detected in the sixties (e.g., Friedman et al. 1967). 
Simultaneous X-ray and radio observations, performed by 
White et al. (1978), did not show the presence of signifi- 
cant radio emission, but revealed correlated X-ray intensity 
and spectral variations similar to those already observed from 
Sco X-1. Observation of the source with the Gas Scintillation 
Proportional Counter (2-20 keV) on-board EXOSAT (White et 
al. 1986) provided an accurate determination of the source con- 
tinuum spectrum, which was fit in the context of the Western 
model with best fit parameters /cTbb ^ 1 keV, radius of the 
BB emitting sphere i?bb ^ 15 km, T ~ 0.6 and Ec 4.5 keV. 
Also an emission line at ^ 6.5 keV with equivalent width (EW) 
of ~110 eV was detected. The data analysis of the same ob- 
servation, but with the Medium Energy Proportional Counter 
Array, showed that the 1-30 keV source spectrum could be 
fit with a BE (fcTbb ~ 1-2 keV) plus the Comptonization 
model proposed by Sunyaev & Titarchuk (1980), with elec- 
tron temperature kT^ ^ 3 keV and optical depth of the elec- 
tron cloud of T 13 (White et al. 1988). The column den- 
sity along the source direction was determined by Vrtilek et 
al. (1991b) with Einstein (0.2-20 keV), who found a value of 
Nu 1.5 X 10^^ cm~^, assuming a simple BB (kTbt ^1.5 
keV) as best fit spectrum. Hoshi & Asaoka (1993) observed the 
source with Ginga in both the HB and the NB. By fitting the 
1-20 keV X-ray spectrum with the first version of the Eastern 
model (BB + DBB), these authors found that the inner disk tem- 
perature and the projected radius of the DBB did not change 
significantly from the NB to the HB, while the BB tempera- 
ture gradually increased from to 2.1 keV to 2.7 keV. DS2000 
analyzed the 0.1-100 keV X-ray spectrum of GX 17H-2 with 
BeppoSAX, discovering a transient hard X-ray tail in addition 
to the persistent continuum. The latter was well described by 
a BE plus COMPTT model, with best fit parameters fcTbb ^ 
0.6 keV, i?bb ~ 40 km, temperature of the seed photons to be 
Comptonized fcT, ^ 1 keV, electron temperature fcTg ^3 keV 
and optical depth of the Comptonizing cloud r ^ 10, assum- 
ing a spherical geometry. Additionally, they found evidence for 
an emission line with energy centroid E\ ^ 6.7 keV and EW 
~ 35 eV, and an absorption edge with £^cdgc ~ 8.6 keV and 
optical depth Tacdge ^ 3 x lO^'^. Type I X-ray bursts from 



GX 17h-2 were observed with the Einstein (Kahn & Grindlay 
1984), HAKUCHO (Tawara et al. 1984) and EXOSAT (Sztajno 
et al. 1986) satellites. 

Despite the X-ray brightness, the optical counterpart of 
GX 17h-2 has not yet been unambiguously identified. The dif- 
ficulties arise mainly because of the position of the source near 
the Galactic center (l — 16°. 4, b = 1°.3), where the chance of 
superposition to unrelated field objects is very high. 

In this paper we report on the results obtained from two 
BeppoSAX observations of GX 17h-2, performed before the 
October 1999 observation reported by DS2000. In Sect.Qwe 
describe the observations and data analysis, in Sect. |3] we 
present the results of the timing and spectral analysis, in Sect.|4] 
we discuss the results and in Sect.|5lwe draw our conclusions. 

2. Observations and data analysis 

We observed the source twice, on 1997 April 3 and 21, with the 
Narrow Field Instruments (NFIs) on board BeppoSAX (Boella 
et al. 1997a). The NFIs include a Low-Energy Concentrator 
Spectrometer (LEGS, 0.1-10 keV; Parmar et al. 1997), three 
Medium-Energy Concentrator Spectrometers (MFCS, 1.5- 
10 keV; Boella et al. 1997b), a High Pressure Gas Scintillation 
Proportional Counter (HPGSPC, 4-120 keV; Manzo et al. 
1997), and a Phoswich Detection System (PDS, 15-300 keV; 
Frontera et al. 1997). Both observations were performed with 
all three MFCS units (MFCS unit 1 failed on 1997 May 6). 
Table ^reports the observations log along with the on-source 
exposure times of each NFI and the mean 1.8-10 keV and 13- 
200 keV count rate measured by all the three MFCS units and 
by the four PDS units, respectively. During all pointings the 
NFIs worked nominally and the source was detected by all of 
them. Good data were selected from intervals when the NFIs 
elevation angle was above the Earth limb by at least 5° and, for 
the LEGS, during the dark earth. The SAXDAS 2.0.0 data anal- 
ysis package was used for processing the LECS, MFCS and 
HPGSPC data, while the PDS data reduction was performed 
using the XAS package v2.1 (Chiappetti & Dal Fiume 1997). 
The LECS and MFCS spectra were extracted from a region of 
4' radius centered on the source position. The background spec- 
trum for the LECS and the MFCS was determined using the 
standard blank field spectra as reported by Fiore et al. (1999), 
while for the HPGSPC and the PDS the background was esti- 
mated with the rocking collimator technique. The spectra of all 
NFIs were rebinned to reduce the oversampling in the energy 
resolution and to have a minimum of 20 counts in each bin, 
in order to reliably use the statistics for spectral fitting. For 
each instrument we performed the spectral and temporal anal- 
ysis in the energy range where the response function is well 
determined; on the basis of this extension,the energy range is 
0.4^.0 keV for LECS, 1.8-10 keV for MFCS and 8-30 keV 
for HPGSPC, while for PDS the range is 13-120 keV in the 
first observation and 13-40 keV in the second one as the source 
was not detected beyond 40 keV. We performed spectral analy- 
sis over four different regions, named I, II, III and IV, according 
to the 1.8-10 keV MFCS count rate being 180 s"\ between 
- 180 and - 185 s^^, between - 185 and - 195 s^^ and <; 
195 s^^ (see Fig.Q. We note that region I is entirely covered 
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Obs. 


Start Time (UT) 


End Time (UT) 


LEGS 


MEGS 


HPGSPG 


PDS 


1.8-10 keV MEGS 


13-200 keV PDS 








ks 


ks 


ks 


ks 


count rate (s^^) 


count rate (s^^) 


1 


1997 Apr 03 00:27:52 


1997 Apr 03 06:00:06 


2.3 


10.7 


4.1 


3.4 


176 


35 


2 


1997 Apr 21 02:29:00 


1997 Apr 21 06:13:55 


1.7 


6.5 


2.4 


2.3 


190 


24 



by the points of the first observation while the other regions 
are covered by the second observation. Only few points of the 
second observation fall in region I, so we could not extract sta- 
tistically significant spectral information in this region during 
the second observation. For the latter we then extracted both 
the average spectrum and the separated spectra for regions II, 
III and IV. 

We used the XSPEC vU.2.0 package to fit the multi- 
instrument energy spectra. In the broad band fits, normaliza- 
tion factors were applied to LEGS, HPGSPG and PDS spec- 
tra following the cross-calibration tests between them and the 
MEGS, constraining them within the allowed ranges (Fiore et 
al. 1999). Interstellar photoelectric absorption, modeled using 
the cross sections implemented in XSPEC and with solar abun- 
dances given by Anders & Ebihara (1982), was included in all 
spectral models used. Uncertainties in the parameters values 
obtained from the spectral fits are single parameter errors at 
90% confidence level, while upper limits are reported at 2a 
level. For the evaluation of the source luminosity we have as- 
sumed a distance of 7.5 kpc (Penninx et al. 1988) 

3. Results 

3.1. Temporal properties 

The source 1.8-10 keV light curve of the two observations 
is shown in Fig. [2 A flux variability on time scales of hours 
is apparent and is confirmed by the fit with a constant func- 
tion, which gives x^/dof = 436/104. Unfortunately, the pres- 
ence of several gaps in the light curve prevents us to investi- 
gate whether the source shows periodic flux modulations. The 
HID of the two observations, obtained defining as intensity the 
1.8-10 keV MEGS count rate and as hard colour the count rate 
ratio in the 6.5-10 keV and 5-6.5 keV energy bands, respec- 
tively, is shown in Fig. |2l As it can be seen, the data points 
of the first observation are all clustered in the left side of the 
diagram, while those of the second observation, reflecting the 
source intensity decrease (from ~ 200 s^^ to ~ 180 s~^, see 
Fig- [D' trace an extended path (regions II, III and IV) that, 
on the left side, smoothly merges with the data points of the 
first observation. However from the shown HID, it is difficult 
to establish which branch of the Z-track our data correspond 
to. To solve this issue, we plot in the lower panel of Fig. |2] 
both our data points and those obtained with BeppoSAX by 
DS2000, when the source clearly traced the HB and NB. As 
DS2000 adopted different energy ranges to produce the HID, 
we used the MEGS event files of their observations from the 
BeppoSAX ^MbYic data archive to derive a new HID in the same 



energy ranges adopted by us. From Fig. |2l it is possible to 
see that GX 17+2, during our observations, traced almost the 
entire HB albeit in a slightly shifted position with respect to 
that observed by DS2000. Secular shifts of the Z pattern in 
the HID of GX 17+2 have been previously observed to oc- 
cur on time scales of several months by RXTE (Homan et al. 
2002) and EXOSAT et al. 1997). On the others hand, 

both RXTE and EXOSAT observations revealed strong stability 
(within ^ 2%) of the hard colour, and this confirms the HB 
nature of the source during our observations. 

We also investigated the source eiTatic time variability in 
the 1.8-10 keV energy band (MEGS data), by deriving the 
power spectral density (PSD) in the 8 x 10^"'-50 Hz frequency 
range, obtained by averaging data intervals of 131072 points 
with a binning time of 10~^ s, for the first and second obser- 
vation respectively. Unfortunately, because of the low collect- 
ing area of the telescopes, the statistical quality of the derived 
PSD did not allow us to resolve the broad band noise com- 
ponents typical of Z sources or to detect QPOs observed with 
other satelhtes (EXOSAT, Hasinger & van der Klis 1990;RXTE, 
Homan et al. 2002). We thus only derived the integrated root 
mean square fractional variation of the 1.8-10 keV flux which 
is 7.1% ± 1.3% in the first observation and 6.3% ± 2.4% in the 
second one. 

3.2. Spectral properties 

Following the results obtained on GX 17+2 by DS2000, we 
first attempted to fit the source spectra with the BB + COMPTT 
continuum model plus a fluorescence Fe line and an absorp- 
tion Fe K-edge. The best fit results are shown in Fig.|3j where 
the time-averaged count rate spectra of the observations 1 and 
2 are shown along with the residuals to this model. As can be 
seen, for the second observation (Fig. O the model describes 
quite well the data, with no apparent systematic residuals, even 
though the is not satisfactory (x^/dof = 192/141). This 
could be due to some spectral evolution of the source with time 
(see Sect. l3.3.2l i. In addition, no significant K-edge is required 
by the data (see Table |2j. Unlike the second observation, the 
above model does not provide a good description of the data 
of the first observation, with a x^/dof = 222/144, mainly 
due to the residuals to the best fit model above 30 keV (see 
Fig. O. It is sufficient to add, as also done by DS2000, a PL 
component to the above spectral model to get a very good fit 
(X^/dof = 143/142). The best fit results are shown in Table 
121 A feature of the data obtained from both observations when 
they are fit with the above model (with or without the PL) is 
that two equally acceptable solutions for the values of the BB 
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Fig. 1. MECS 1.8-10 keV light curves of the first (upper panel) and second {lower panel) observation of GX 17+2. The bin size 
in both cases is 100 s. Time is given in seconds from the beginning of each observation (see TableQ. 



temperature kTbb and temperature of the seed photons kTg are 
found (see Table |2j: one with fcTbb < kTg and the other with 
kTbb > kTs, while the other parameter values remain substan- 
tially unchanged within errors. 

Finally, to test the physical hypothesis of the three-layered 
atmospheric structure in accretion disks suggested by Zhang et 
al. (2000) for the black hole candidates (BHC) GRO J 1655-40 
and GRS 1915+105, in the first observation we replaced the PL 
component with a further COMPTT model with higher electron 
temperature, which we label as kT^. Approximating the seed 
photons temperature of this hotter electrons to the temperature 
of the electrons of the underlying warm layer (an assumption 
justified by the fact that the Compton parameter y of the warm 
electrons is '--^ 3, so that their emerging spectrum is not far from 
a Wien-like), namely keeping kT^^ — kTc, does not give any 
stable solution. A good fit is obtained by keeping free all the 
model parameters (x^/dof = 149/140), but in this case, unlike 
the results by Zhang et al. (2000), the seed photons tempera- 
ture of the hotter COMPTT component (kT^ ~ 1.3 keV, with a 
corresponding kT^ ~ 75 keV and ^ 0.05) is quite different 
from the electron temperature 8 keV) of the warm electrons. 

3.3. A physical model for the hard X-ray tail 

In order to investigate the physical origin of the hard X-ray 
tail detected during the first observation of GX 17+2, we used 
the Comptonization model worked out by Coppi (1999, EQ- 



PAIR in XSPEC), which has the advantage of being general 
and can be used to fit the source spectra of both observa- 
tions. In this model, non-thermal electrons with energy spec- 
trum oc (where 7 is the Lorentz factor) are injected in 
the plasma cloud at a rate (x 7^^. These electrons lose energy 
because of Compton, Coulomb and bremsstrahlung processes 
and thus establish a steady-state distribution. At high energies 
(for 7 > 7inin up to a certain value of 7max) the distribution 
is non-thermal (power-law like), but at low energies it joins 
the Maxwellian distribution and the final result of this process 
is a hybrid thermal plus non-thermal plasma. The temperature 
of the thermal distribution, kT^, is self-consistently computed 
by means of energy balance and is not a free parameter of the 
model. The plasma properties are characterized by the powers 
L supplied to the different components. A compactness param- 
eter is introduced and defined as: 



where R is the characteristic size of the region and ctt the 
Thomson cross-section for scattering. The compactness of the 
input soft photons is parametrized by t^. The power suppUed 
to the Comptonizing electrons is parametrized as £h=^th+ ^nth^ 
where £th is the power supplied to heat thermal particles, while 
Inth is the power going into accelerating non-thermal particles. 
If ^nth— 0, then one has a purely thermal (Maxwellian) distribu- 
tion. EQPAIR can also test the presence of Compton reflection 
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Fig. 2. Left panel: HID of the two BeppoSAX observations of GX 17+2. Diamonds: first observation. Open squares: second 
observation. The four count rate regions are also shown. Each point corresponds to a bin size of 100 s. Right panel: source HID 
of both our observations (same symbols as in the left panel) and that performed by DS2000 {crosses). The bin size of the 1997 
October data points is 400 s. A typical error bar is shown on the top right of each panel. 
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Fig. 3. Count rate spectra of the first observation {left panel) and of the second observation {right panel) of GX 17+2 and best 
fit folded model WABS(BB + COMPTT + GAUSSIAN) with fcTbb < kT^. In the first observation, an absorption edge was also 
included. Below each panel are shown the residuals to the model in units of a. The last two PDS points of the second observation 
have S/N < 3, so they were not included in the fit and are reported as 2a upper limits. Similar results are obtained for the same 
model but with fcTbb > kT^. 



of the primary photon spectrum from a neutral or a ionized 
medium (Magdziarz & Zdziarski 1995) allowing for different 
metal abundances with respect to the solar ones. 

3.3.1 . First observation 

We first assumed that all the power goes in heating the thermal 
particles (^nth/^h= 0, purely thermal Comptonization) with a 
simple BE spectrum of the seed photons. We added to EQPAIR 



also a direct BE component plus a Fe K fluorescence line in the 
source model, as discussed in Sect. 13.21 This model well de- 
scribes the data below 30 keV but, as in the case of COMPTT, a 
systematic excess is again observed in the residuals above this 
threshold (see Fig.l^. We thus left free to vary the ^nth param- 
eter: since the statistics of the data do not allow us to simulta- 
neously constrain 7niin, 7max and p in the fits, 7inax was frozen 
at value of 12 (the is poorly sensitive to 7max in the range 
from 2 to 50). The chosen value of 7niax takes into account the 
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Fig. 4. Residuals in units of a to the model WABS(BB + EQPAIR + GAUSSIAN), assuming purely thermal Comptonization, for the 
first {left panel) and for second (right panel) observation of GX 17+2. 

Table 2. Best fit parameters of the multi-component model WABS(BB + COMPTT + GAUSSIAN + K-edge + PL) for the first 
and second observation of GX 17+2. The two solutions of the model are reported. The F-test gives the probability of chance 
improvement of the for the addition of the PL. The ratios ibb/^tot and Lpi/Ltot are computed in the 0.4-200 keV energy 
range. 





Obs. 1 


Obs. 2 


Obs. 1 


Obs. 2 




I 


II+III+IV 


1 


11+111+lV 


Parameter 


Solution 1 


Solution 2 


kTib (keV) 
Rbb (km) 
kTs (keV) 
kTe (keV) 

T 

r 

F-test 


9 c;+0.2 

3511 

1-^^-0.08 

3.4lg:i 

9 O+O.l 
9 ,+1.2 

3 X 10-1* 


9 9+0.1 
^■^-0.2 

0.6ltH^3 
40tl 

r.ua_o.o4 

T Q+O.l 

-J-J-0.1 

10 5+0-4 
[2.8] 

f, 7+0.6 
^•'-0.6 

0.11 


9 c+0.2 

1 7Q+0.10 

J^- '^-0.16 
T 7+0.6 
-^•'-0.6 
r, ;ro + 0.02 
^J-JO-0.02 

3.5lS]:i 

12.0t°-j 

T O+O.l 

9 0+1.2 
^•■5-1.3 

5 X 10-12 


9 q+0.1 
^■"J-0.1 

1 qO + 0.05 

i.OO_o.o2 

7 4+0.5 

U.00_o.o3 

q 1 c+o.oe 

•J-^>J-0.04 

13.0«J 
[2.8] 

^■^-0.6 

0.03 


El (keV) 
(Ti (keV) 

EWi (eV) 
Sedge (keV) 

Tedge(10-2) 


f, 7+0.1 

i+°-i 

S T+0.4 
°-^-0.3 

^•-'-0.7 


f: 7+0.1 
"■'-0.1 

[8.2] 
< 2.0 


f, 7+0.1 

-0.1 

i+° i 

q + 1.2 

27t« 

n O+0.3 

0-^-0.2 

3.3lli 


6.7l°o:i 
o.2i°:l 

4 ts+i-'5 

^■0-1. 

31±^i 
[8.3] 
< 1.8 


Lhh/ Ltot 

Lpi/ Ltot 

^ 0.4-30 kcV 
^ 30-200 kcV 

X^/dof 


0.16 
0.30 
1.87 
0.02 
143/142 


0.17 
0.13 
1.62 
0.01 
188/140 


0.10 
0.32 
1.91 
0.02 
140/142 


0.15 
0.18 
1.69 
0.01 
189/140 



^ In units of 10^^ cm-^. 

^ Photons keV-i cm-^ s-i at 1 keV. 
Total photons in the line in units of IQ-'^ cm-^ s 
Unabsorbed luminosity in units of 10"^** erg s-i. 



energy range in which the excess over thermal Comptonization 
is located and the temperature of the seed photons found from 
the fit. As can be seen from Table|3] the agreement between the 
model and the data is excellent. The probability given by the F- 
test that the reduction obtained from the inclusion of a non- 
thermal component in the Comptonizing cloud is due to chance 
is very low (see Table|3}. The presence of a Compton reflection 
component, assuming solar abundances and a low inclination 
angle i = 30° (GX 17+2 is a Sco-like source, Kuulkers & van 
der Klis 1996), is only marginally detected (Ax^'^ 7, probabil- 
ity of chance improvement of 3%) with best fit value of the 
covering factor n/2TT ^ 0.04. The ionization parameter value 



we find is ^ 800 erg cm s"^. An important issue to point out 
is that, unlike COMPTT, with EQPAIR there is not a double so- 
lution as concerns the direct BB and seed photons temperature 
values. 



Replacing the direct BB component with a DBB, we still ob- 
tain an acceptable, albeit a bit worse, result (x^/dof = 156/142). 
If, on the other hand, we assume a DBB seed photon spectrum 
in EQPAIR plus a direct BB model, the fit is even worse (x^/dof 
= 163/142). 
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Table 3. Best fit parameters of the multi-component model wabs(bb + eqpair + Gaussian) for the first and second observation 
of GX 17h-2. The F-test gives the probability of chance improvement of the switching from thermal to hybrid Comptonization. 
The compactness parameters is, £h and ^„th are defined in Sect. 13. 31 The ratio ibb/^tot is computed in the energy range 0.4-200 
keV. 



Parameter 


Obs. 1 

T 
1 


Obs. 2 
TT-i-TTT-i-TV 

llTlilTl V 


iVg 
kUb (keV) 
Rbb (km) 
kTs (keV) 
kn (keV) 

is 

(■nth/t-h 
T 
P 
7min 

F-test 


o.581°:E;? 

33+^ 

3.5 

7600+^"°°° 

r, 87+0.02 
"•°'-0.02 

^- A ^ _0 02 

4+1.5 

1 7+1-^ 

1- '-0.6 

i-i-'-o.is 
[12] 

4 X 10~" 


, 07+0.07 

3111 

0.95lH'5 
3.5 
[7600] 

< 0.04 
04+0.1 

[1.7] 
[1.15] 
[12] 
0.07 


El (keV) 
o-i (keV) 

EWi (eV) 


-,7+0.1 

o.2«:} 

i2tl 


-0.1 
2+"-^ 

4-o_o.g 
33115 


Lhb/Ltot 
^ 0.4-30 keV 
^ 30-200 kcV 

X^/dof 


0.11 
1.44 
0.019 
143/142 


0.11 
1.46 
0.006 
185/141 


In units of 10 


22 -2 

cm . 





^ Total photons in the line in units of 10 cm ^ s 
Unabsorbed luminosity in units of 10'^* erg s^^. 



Table 4. Same as in Table|3]but for regions II, III and IV of the second observation as reported in Fig.|2]and assuming thermal 
Comptonization. 



Parameter 


11 


Ill 


IV 


iVg 
kTbb (keV) 
Rbb (km) 
kTs (keV) 


2 03+0.11 

^■^J-0.05 

61+° '^i 

371^3 
, ,7+0.02 

'-0.07 


2 0^5+0.10 

^■^J-0.13 
5^+0 09 
"•-'■-'-0.06 
^7 + 10 
■^^-10 
n o;r+0.05 
U.63_o.08 


1 03+0.18 
l^"--'-0.14 

0.55lH«6 
291^ 

ngn+O.lO 
'-'•^"-'-0.04 


kTc (keV) 


4.0 


3.5 


3.1 


T 


[7600] 
47+0 01 

"•^'-0.03 
f: 7+0.2 
''•'-0.5 


[7600] 

^■(^^ 


[7600] 
50+0 02 

"•J7_0.02 

8.9l°J 


El (keV) 
m (keV) 

EWi (eV) 


f. ^+01 

o.2i»:^ 

5.01^1 

3^+25 
■^■^-16 


f, -7+0.2 
"•'-0.2 

0.3l°i 

6.5lt.? 
44+33 


6^8lS:l 

i+°-2 

'■'•1-0.1 

I- 0+2.4 
^••J-1.8 

3f,+ 16 


Lhh/ Ltot 


0.17 


0.07 


0.07 


^ 0.4-30 kcV 


1.47 


1.50 


1.48 


^ 30-200 kcV 


0.01 


0.007 


0.004 


X^/dof 


145/142 


133/142 


145/142 


" In units of 10' 


22 cm-2. 






Total photons in the line in units of 10 ^ cm ^s ^. 




Unabsorbed luminosity in units of 10^^ erg s"^. 
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Fig. 5. Unabsorbed EF{E) spectrum of the first {left panel) and of the second {right panel) observation of GX 17+2 along with 
the best fit model WABS(BB + eqpair + Gaussian). A hybrid spectrum of the Comptonizing electrons is assumed. The 2a 
upper limit of the last two PDS points of the second observation is also reported. Below each panel, the residuals between the 
data and the model in units of a are shown. Different line styles represent the single components. Dotted: BE. Dotted-dashed: 
EQPAIR. Dashed: Gaussian emission line. 



3.3.2. Second observation 

As reported in Sect. 13.21 during the second observation there 
was slight evidence of a spectral evolution of the source, sug- 
gested by the somewhat worse fit of the COMPTT model to 
the data, in spite of the fact that no systematic residuals to 
the best fit model are apparent (see Fig.|3}- This result is con- 
firmed also by using the EQPAIR model, assuming purely ther- 
mal Comptonization, to fit the whole observation (x^/dof = 
186/142). In this case we also find that the soft compactness 
Is is unconstrained so, during the fit, we fixed it at the best fit 
value of the first observation. In order to put an upper limit on 
the power supplied in accelerating non-thermal electrons, we 
allowed the ^nth/^h parameter to vary, by keeping 7,jii,i and p 
frozen at the best fit values of the first observation. The results 
are reported in Table|3] As it can be seen, the power going in ac- 
celerating non-thermal electrons is consistent with zero and its 
2(7 upper limit is lower than the value of the first observation. 
It is worth pointing out that this upper limit value on IntiJ^h 
is obtained even including in the fit the last two PDS points 
(see Fig. |3j,in order to have the same energy coverage of the 
first observation. The 2a upper limit on the 40-120 keV source 
flux, assuming as best fit model the one reported in Table |3] is 
~ 7 X 10~^^ erg cm~^ s~^ (see Fig.|5} a value below that found 
in the first observation (^ 1.0 x lO^^*' erg cm^^ s^^). We sub- 
sequently performed spectral analysis separately for the three 
regions of the HID labeled in Fig.|2j the results are reported in 
TableS 

4. Discussion 

Some relevant results have been derived from our observa- 
tions of GX 17+2 with BeppoSAX. The persistent continuum 
X-ray emission of GX 17+2 can be described by several two- 
component models. When we adopt the comptt model plus 



a BB, we find that two solutions are possible, which are remi- 
niscent of the well known Eastern (Mitsuda et al. 1984, 1989) 
and Western (White et al. 1986, 1988) models. Indeed, looking 
at the results of Tabled for solution 1, we can interpret the ^ 
0.6 keV BB as the direct emission from the accretion disk, with 
Comptonization mainly occurring in the hotter boundary layer 
(fcTg ^ 1.2 keV, Eastern Model scenario). But the same spec- 
trum can be described (solution 2) by a ~ 1.5 keV BB plus a 
Comptonization region of cooler photons (fcTg ^ 0.6 keV). In 
this case the direct BB emission can be interpreted as coming 
from the boundary layer, while Comptonization takes place in 
the disk (Western model scenario). From a statistical point of 
view, we do not find strong reasons for favouring one solution 
over the other; however, the former solution is the only one we 
find when COMPTT is replaced by EQPAIR. 

The soft compactness estimated by EQPAIR (see Table |3} 
is much higher than that reported in the case of BHCs (e.g., 
Frontera et al. 2001; Gierliiiski & Done 2003); this is how- 
ever not surprising given that in neutron star systems, unlike 
BHCs, the source environment in the inner disk region has 
an additional source of soft luminosity coming from the neu- 
tron star. Apart from this crucial difference, related to the na- 
ture of the compact object, there are some physical character- 
istics which link in some way GX 17+2 to BHCs. First, the 
spectral index of the non-thermal injected electrons {p ^ 1 .7) 
is consistent, within errors, with that found, e.g., in Cyg X- 
1 (Frontera et al. 2001) and GRS 1915+105 (Zdziarski et al. 
2001). Moreover, the spectral shape of the source in the EF{E) 
diagram (see Fig. |5} strongly resembles that of BHCs in the 
so-called ultrasoft state, when Lx ^ -^Edd (e-g-, GX 339-4, 
GRS 1915+105, XTE J 1550-564, see the review by Zdziarski 
& Gierliriski 2004), where the peak of the energy emission oc- 
curs below 10 keV. The similarity with GRS 1915+105, in par- 
ticular, results even more intriguing looking at the X-ray spec- 
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tral parameters found by Zdziarski et al. (2001) by fitting the 
source RXTE/OSSE 3-600 keV X-ray spectrum with EQPAIR: 
the authors indeed found kT^ --3.6 keV, r ~ 4.4, kT^ 
1.3 keV and ^nth/4h^ 0.1. Only the ratio ^h/4was lower 
0.3) than that found in GX \l+2 by us. As Z sources stably 
accrete at near-Eddington rate, these similarities lead to sug- 
gest that Comptonization/acceleration processes in X-ray bi- 
nary systems in this phase do not strongly depend on the pres- 
ence or not of a solid surface at the inner boundary of the ac- 
cretion disk. 

An important result of our study concerns the transient high 
energy X-ray tail. We confirm its presence when the source is 
in the upper left HB of the Z diagram, as already observed by 
DS2000. But, unUke DS2000, who found that the hard X-ray 
excess gradually decreased as the source moved towards right 
in the HB, completely disappearing when the source achieved 
the NB, we find that this disappearance occurs when the source 
is still in the HB, in particular during its transition from re- 
gion I to region II of the HID (Fig. |2li. The hard X-ray ex- 
cess is described by a PL with the same photon index found by 
DS2000 (r ~ 2.7, see Table|2ji. In order to gain some insight on 
the physical origin of the hard X-ray component, we replaced 
the COMPTT model with the Comptonization model EQPAIR of 
Coppi (1999). In this model, in addition to a thermal popula- 
tion of electrons, a non-thermal electron tail can be included. 
Assuming a purely Maxwellian distribution of the electron en- 
ergies, the hard X-ray excess is still apparent above the BB plus 
EQPAIR model (see Fig.l^. Only assuming a thermal plus non- 
thermal plasma the EQPAIR model gives a good description of 
the high energy excess. To our knowledge, this is the first time 
that the transient hard X-ray tail of a Z source (see Sect. ^ is 
described in terms of hybrid Comptonization. 

An important question is whether the two electron pop- 
ulations (thermal and non-thermal) are physically connected. 
Unfortunately the present data do not allow to answer this ques- 
tion. We consider two cases separately. If the two electron pop- 
ulations are not physically related, a likely origin of the non- 
thermal electrons is a relativistic jet escaping the system. In re- 
cent years, relativistic jets have revealed to be a quite common 
property of LMXBs, providing a link between these systems 
and AGNs (for a review see, e.g.. Fender 2002). In this ther- 
mal plus non-thermal scenario, the Comptonized spectrum is 
produced by a thermal corona (likely located above the inner 
disk close to the boundary layer) plus a relativistic non-thermal 
electron jet, with the soft seed photons mainly coming from the 
inner disk. The energy of the seed photons (^ 1 keV, see Table 
O is indeed consistent with the value expected for the inner 
disk temperature of LMXBs hosting a neutron star (see, e.g., 
the review by Barret 2001). Thus we could tentatively identify 
the decrease of the power supplied to non-thermal electrons 
(^nth/^h drops from 0. 1 to less than ^ 0.04 from the first to 
the second observation), as a switching-off of the mechanism 
producing the jet. 

On the other hand, if the two electron populations are physi- 
cally related, a possible origin of the non-thermal component is 
in processes of magnetohydrodynamic turbulence occurring in 
the corona (e.g., magnetic reconnection), in analogy with what 
observed in the solar corona (Crosby et al. 1998, see also the 



review by Coppi 1999). In this case particle acceleration cannot 
be 100% efficient in the production of non-thermal high energy 
particles (e.g., Zdziarski et al. 1993). This can occur only above 
a certain momentum threshold when the acceleration time scale 
for a particle becomes significantly shorter than its correspond- 
ing thermalization (or energy-exchange) time scale. Below this 
threshold particles are tightly coupled and their acceleration 
goes mainly in a bulk heating of the thermal plasma compo- 
nent. 

We have also explored the physical scenario proposed by 
Zhang et al. (2000) for GRO J1655-40 and GRS 1915H-105, in 
which a much hotter, optically thin thermal corona is located 
above a warm layer with lower temperature and high opacity. 
This scenario seems not suitable for GX 17h-2. In fact, in this 
hypothesis, one would expect that the seed photons of the hot- 
ter layer are those coming from the underlying one, and thus 
should have a temperature similar to that of the colder electrons 
{kT^ - 3-5 keV). This is indeed what Zhang et al. (2000) 
found with ASCA/RXTE data in the case of GRO J 1 655-40 and 
GRS 1915H-105. However, in our case kT^} significantly differs 
from fcT"' 1.3 keV and ^ 8 keV, respectively), thus ruling 
out the three-layered accretion hypothesis. 

In addition to the origin of the hard X-ray emission, its 
correlation with the position of Z sources in the CD/HID 
is still unclear. As discussed in Sect. ^ the results obtained 
from GX 17H-2, Cyg X-2 and GX 5-1 (but not, however, from 
Sco X-1), could suggest the possibility of an anti-correlation 
between presence of hard X-ray emission and accretion rate 
M in Z sources, in a similar fashion to what observed both in 
atoll sources and BHCs, where the hard spectra are usually ob- 
served during the low inferred M states (e.g.. Barret & Olive 
2002; McClintock & Remillard 2004). Our results on GX 17h-2 
do not contrast with this possibility but, when compared with 
those found by DS2000, they point to the presence of a thresh- 
old, that we thus tentatively identify with a threshold in M, 
above which the hard tail disappears. A definite answer to this 
suggestion could come from very sensitive monitoring of the 
source at high energies in region I of Fig. |2l during the second 
observation, but this was not possible because of the too few 
points falling in that region. 

An important issue is what we mean by M: it appears ev- 
ident that this M is not the one that drives the X-ray luminos- 
ity, since the latter does not change along the track (see Tables 
|3land|3. Thus, if some M increases along the Z-track, it is 
not correlated with the X-ray luminosity: an intriguing way 
out to this problem was proposed by van der Klis (2001) to 
explain the "parallel tracks" phenomenon. He suggested that 
what does change along the track is not the prompt disk ac- 
cretion rate Md{t) but rather the ratio ?/ = Md(t)/ {Md{t)), 
where (il/d(t)) is the long-term A/d(i) average. The luminos- 
ity is supposed to be related to the sum of Md{t) plus (Afd(t)) 
(Lx oc Md(i) + oi{Md{t))) while the curve length S across 
the track is a measure of ij. The threshold above which we see 
the disappearance of the hard tail could thus be a threshold in 
Md(t)/(A/d(i)). It is quite difficult to estabUsh whether this 
hypothesis is correct. In the positive, one has to explain what 
produces the hard tail: the physical process should be strongly 
dependent on what happens in the disk. When the hard X-ray 
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component is no longer detected (regions II, III and IV of the 
HID), we observe a slight decrease of the temperature of the 
scattering cloud, as the source moves from the left to the right 
in the HB. In the van der Klis (2001) model, this is related to 
the variation of M: a higher prompt disk accretion rate is ex- 
pected to push the inner radius towards the neutron star with a 
consequent larger production of seed and hotter photons for the 
corona, which is thus more efficiently Compton-cooled (kTc 
decreases from regions II to IV, see Table 0}. We confirm the 
presence of the Fe K emission line at ~ 6.7 keV, while the 
edge detection (first reported by DS2000) is not so evident and 
should be treated carefully. Indeed we find that the EQPAIR 
model well fits the data (see Table |3j and neither Compton- 
reflection nor absorption edge are required, showing that the 
presence of such components in the spectra are very sensitive to 
the continuum model adopted. One should be careful in claim- 
ing the presence of these components in the case of GX 17h-2. 

5. Conclusions 

We have investigated two BeppoSAX observations of the Z 
source GX 17h-2 performed on 1997 April 3 and 21. The per- 
sistent continuum X-ray emission can be described by two- 
component models, which support the well known accre- 
tion scenarios suggested by the Western and Eastern model. 
However, in the first observation, when the source was in the 
upper left HB of the HID diagram, we observed a hard X-ray 
tail (up to 120 keV), which disappeared in the second obser- 
vation performed two weeks later, when the source moved to 
the right in the HB. Our results confirm those obtained by 
DS2000, but provide a stronger evidence for the presence of 
a threshold in the accretion rate M in the HB, above which 
the intensity of the hard tail is below the instrument sensitiv- 
ity or disappears. Given that the X-ray luminosity is not corre- 
lated with the position of the source in the HID, M is not the 
one we ascribe to X-ray luminosity but has a somewhat dif- 
ferent meaning. We propose as working hypothesis the model 
by van der Klis (2001), in which the X-ray luminosity is pro- 
portional to the prompt disk accretion rate plus its own long- 
term average, while motion along the Z track depends on the 
ratio of these two quantities. To testify whether this hypothe- 
sis is correct, more observational and theoretical work is re- 
quired. The hard X-ray excess observed in the first observa- 
tion can be described by Comptonization of soft seed photons 
(kTs ^ 1 keV) off non-thermal electrons injected with spec- 
tral index p 1.7. In addition, a thermal component described 
by a BB with fcTbb ^ 0.6 keV is observed. An emission line 
at ~ 6.7 keV is found in both observations. The origin of the 
non-thermal electrons cannot be established on the basis of the 
BeppoSAX data alone. We have two equally acceptable options: 
electrons from jets or electrons produced in magnetohydrody- 
namic turbulence processes occurring in the corona. Only the 
determination of the source spectrum at higher energies (<; 200 
keV) can provide more constraints on the origin of the high 
energy excess. A step forward could be done by INTEGRAL 
through a long monitoring observation of the source. Also high 
energy resolution observations below 10 keV, performed with 
Chandra or XMM-Newton, could be very helpful in clarify- 



ing the origin of the Fe emission line and to determine unam- 
biguously the presence/absence of an absorption edge. On the 
other hand, simultaneous X-ray/radio observations would test 
the possibility of a correlation between hard X-ray and radio 
emission in this source, with presence or not of relativistic ex- 
panding lobes. 
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